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Abstract. We present observations with the adaptive op- 
tics assisted integral field spectrograph OASIS of the M 31 
double nucleus in the spectral domain around the Calcium 
triplet at a spatial resolution better than 0'.'5 FWHM. 
These data are used to derive the two-dimensional stel- 
lar kinematics within the central 2". Archival WFPC2/HST 
images in the F300W, F555W and F814W bands are re- 
visited to perform a photometric decomposition of the 
nuclear region. We also present STIS/HST kinematics ob- 
tained from the archive. 

The luminosity distribution of the central region is well 
separated into the respective contributions of the bulge, 
the nucleus including PI and P2, and the so-called UV 
peak. We then show, using the OASIS kincmatical maps, 
that the axis joining PI and P2, the two local surface 
brightness maxima, does not coincide with the kinematic 
major-axis, which is also the major-axis of the nuclear 
isophotes (excluding PI). We also confirm that the veloc- 
ity dispersion peak is offset by ~ 0'.'2 from the UV peak, 
assumed to mark the location of the supermassive black 
hole. The newly reduced STIS/HST velocity and disper- 
sion profiles are then compared to OASIS and other pub- 
lished kinematics. We find significant offsets with previ- 
ously published data. Simple parametric models are then 
built to successfully reconcile all the available kinematics. 

Wc finally interpret the observations using new N-body 
simulations. The nearly keplerian nuclear disk of M31 is 
subject to a natural m = 1 mode, with a very slow pat- 
tern speed (3 km/s/pc for Mbh = 710 7 Mq), that can 
be maintained during more than a thousand dynamical 
times. The resulting morphology and kinematics of the 
mode can reproduce the M 31 nuclear-disk photometry 
and mean stellar velocity, including the observed asym- 
metries. It requires a central mass concentration and a 



cold disk system representing between 20 and 40% of its 
mass. Such a slow mode could be excited when interstel- 
lar clouds from the more external gaseous disk infall to- 
wards the centre. Nuclear disks formed from accreted gas 
are possible candidates for the precursors of these types 
of structure, and may be common in central regions of 
galaxies. 

Key words: Galaxies: individuals: M 31; Galaxies: kine- 
matics and dynamics; Galaxies: nuclei; Galaxies: photom- 
etry; Instabilities 



1. Introduction 

The role of supermassive black holes (hereafter SBHs) as 
power sources for active galactic nuclei is now widely ac- 
cepted. The number of SBH candidates is growing rapidly 
and with more than 20 objects with good black hole mass 
estimates, one can start to search for statistical relation- 
ships between the black hole mass and the surrounding 
galaxy (de Zeeuw 2000] and references therein). Prelimi- 
nary studies indicate that the SBH mass is correlated with 



the luminosity of the host galaxy (Magorrian et al. 1998) 
and, with much less scatter, with its velocity dispersion 
(Ferrarese & Merritt |2000| , Gebhardt et al. [2000 ). This 



suggests that the SBH mass is of the order of 10 d times 



* Based on observations collected at the Canada-France- 
Hawaii Telescope, operated by the National Research Council 
of Canada, the Centre National de la Recherche Scientifique of 
France, and the University of Hawaii 



the mass of the bulges (Merritt & Ferrarese 2000). If this 
relationship is confirmed, it would have important impli- 
cations on our understanding of the formation of SBHs in 
galaxies. However the statistics is still biased to high SBH 
mass because of observational as well as modeling limita- 
tions. Extending the study to a larger number of objects 
and to lower SBH mass will require a large observational 
effort. 

Detailed study of individual objects is fundamental not 
only to better constrain the SBH mass, but also to pro- 
vide observational support on the origin of a possible re- 
lation between the SBH and global properties of the host 
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galaxy. Theoretical work suggests that SBHs may influ- 
ence the central morphology of galaxies through scatter- 
ing of centrophilic stars (see e.g. the review by Merritt 
1999). Nuclear bars or spirals have been observed in the 



central region of galaxies, and thus proposed as a trigger- 
ing mechanism for the growth of SBHs. Although this does 
not solve the problem of driving the dissipative component 
within the central parsec, it strongly suggests that the hy- 
pothesis of axisymmetry and steady-state widely used to 
model galaxy cores and constrain the SBH mass may not 
be adequate. 

In this paper we will focus our attention on the nu- 
cleus of M 31, which is an ideal case for such a detailed 
study: it is well resolved with a diameter of approximately 
4" (^15 pc), even at ground-based spatial resolution, it 
is bright (my ~ 13), does not suffer from dust absorption 
(contrarily to the Galactic centre) and has a favourable 
inclination (i ~ 70°). There is a long history of kine- 
matical observations of this nucleus, with first evidence 



for a ~ 10 7 M centr al SB H by Kormendy ( ]1988| ) and 



Dressier & Richstone ( |1988| ). It stood out then as one 
of the most convincing cases, since it was inferred from 
stellar kinematics obtained from high signal-to-noise ratio 
spectroscopic observation at 4 pc (~ 1") resolution. How- 
ever the modeling was based on the unrealistic assump- 
tion of spherical symmetry. Upgrading the modeling to an 
axisymmetric (or triaxial) geometry would have been pos- 
sible, but it was already known since the Stratoscope II 



observations (Light et al. 1974), that the light distribution 
was not even symmetric. This asymmetry was resolved in 



an intriguing double structure by Lauer et al. (1993) using 



the HST (pre-COSTAR) WFPC imaging capabilities. In th e 
HST original and post-COSTAR images (Lauer et al. |l998|) , 
the double nucleus appears as a bright peak (PI) offset by 
~ 0'.'5 from a secondary fainter peak (P2), nearly coincid- 
ing with the bulge photometric centre, and the suggested 
location of the SBH inferred from the previous spectro- 
scopic observations (see the discussion in Kormendy & 
Bender 199S). HS T images fro m the far-UV to near-IR 
(King et al. 1995 ; Davidge etal. 1997 ) and long-slit spec- 
tra (Kormendy & Bender 1999] ) all demonstrate that PI 
has the same stellar population as the rest of the nucleus, 
and that a nearly point-like source produces a UV excess 



close to P2 (King et al. |1995J). 

The first available two-dimensional maps of the stel- 
lar kinematics of M 31, obtained by Bacon et al. (1994) 



with the TIGER integral field spectrograph, showed an- 
other unexpected feature: while the stellar velocity field 
is roughly centred on P2, the peak in the velocity disper- 
sion map is offset by ~ 0'.'7 on the anti-Pi side. Further 
HST spectrographic observations were conducted with FOS, 
as well as FOC (Statler et al. 1999| ). In the following we 
will solely refer to the long-slit FOC observations of Statler 
et al. ( 1999] ) as they are the only published data. At the 
HST resolution, the velocity curve presents a strong gradi- 
ent and the zero velocity point is offset by ~ 0'.'16 from P2 



towards PI. The velocity dispersion peak reaches a value 
of a with 440 ± 70 kms -1 (to be compared with the best 
corresponding ground-based value of ~ 248 kms , Kor- 
mendy & Bender |1999[ ). In the Statler et al. ( |1999| ) data, 
the dispersion peak is nearly centred on P2 (within 0'.'06). 
The FOC measurements however suffer from a low signal- 
to-noise ratio (S/N hereafter) (~ 14 per pixel at P1/P2), 
and were obtained via a complex data reduction proce- 
dure to palliate some calibration problems. The shapes of 
the velocity and velocity dispersion curves should thus be 
confirmed with better S/N data . Such data have been ob- 
tained by Kormendy k, Bender (|l999| ) using the SIS/CFHT 
spectrograph although with a significantly lower spatial 
resolution of 0'.'64 FWHM. They confirmed that the dis- 
persion peak is offset by ~ 0'.'2 (in the direction opposite 
to PI) from their assumed velocity centre, and that the 
nucleus is significantly colder than the surrounding bulge. 

On the modeling and theoretical fronts, the situation 
remains uncertain. There is a consensus regarding the ex- 
istence of a SBH of a few 10 7 Mq although published 
values range from 3 to 10 10 7 Mq. Such a black hole mass 
is needed to explain the high value of the velocity disper- 
sion and the fast rotating stellar disk. Note however that 
the precise SBH mass cannot be accurately derived since 
we do not yet have a self-consistent model that can ac- 
count for the observed properties. The SBH is assumed 
to coincide with the centre of the UV peak, near P2, and 
possibly with the recently uncovered hard X-ray emission 
detected by Chandra (Garcia et al. 2000 ). In that scheme 
the position of the velocity dispersion peak should coin- 
cide with the location of the SBH, roughly consis tent with 
the kinematics presented by Statler et al. ( |1999| ), but not 



with the ground-based observations of Bacon et al. (1994) 
and Kormendy & Bender (1999). 

The nature of PI and the observed dynamical pecular- 
ities remain a puzzle. While various possibilities have been 
discussed, only two hypotheses have been studied in some 
detail. Tremaine ( 1995| ) first proposed a model where an 
eccentric disk of stars orbiting the SBH at P2 produced 
the observed accumulation of light at PI. This parametric 
model was able to reproduce the light distribution and 
some of the kinematical features. Emsellem & Combes 



( 1997 ) built the first (and still only) self-consistent mod- 
els for the nucleus of M 31, via N-body simulations. They 
performed simulations of a stellar cluster captured in the 
SBH potential and showed that the available observa- 
tional properties could be reproduced well with properly 
tuned orbital elements for the falling cluster. However the 
timescale for the disruption of the cluster is quite small 
(~ 10 5 years), significantly weakening this scenario. More 
problematic is the observed homogeneity in the colours of 
P2 and PI whose explanation would require some more 
fine tuning of the stellar population of the cluster. 

As emp hasize d by Statler et al. (|1999 ) and Kormendy 
& Bender ( |l999| ), the Tremaine (|1995|) model then re- 
mained the only attractive solution: it naturally predicts 
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the same stellar population for PI and the rest of the nu- 
cleus and seems to be, at least qualitatively, compatible 



with the observed kinematics. Statler (1999) noted that it 
may also help to explain a wiggle observed in the FOC ve- 
locity profile near the location of PI. In Tremaine's origi- 
nal paper, some preliminary suggestions were given for the 
possible formation of such an eccentric disk. These ideas 
have yet to be confirmed by more realistic self-consistent 
models. Such an investigation is also important to deter- 
mine if eccentric disks are common in SBH environments. 

The goal of this paper is twofold: (i) add more obser- 
vational constraints to the photometry and kinematics of 
M 31 nucleus and eliminate the uncertainty regarding the 
location of the velocity dispersion peak; (ii) investigate 
Tremaine's model or alternative solutions in more detail. 

To achieve these goals, we have performed new ob- 
servations with the adaptive optics assisted integral field 
spectrograph OASIS (Section ||) . Although these new kine- 
matical maps have a factor two better spatial resolution 
than previous two-dimensional spectroscopic data (Bacon 
et al. 1994), they cannot compete with the HST spatial res- 
olution. The relative centring and orientation of the vari- 
ous key features observed in the nucleus are also critical. 
We therefore performed a detailed analysis of the photom- 
etry using archival HST/WFPC2 images (Sect. ||). We thus 
combined them with new velocity and velocity dispersion 
profiles derived from the recently released STIS/HST data 
(Section [|). Together, these data provide new constraints 
on the dynamics of M 31 (Sect. ||). We finally interpret the 
observations in terms of slow m = 1 modes and present 
new N-body simulations to support this scenario in Sec- 
tion ||. Conclusions are given in Sect. 0. 

2. HST/WFPC2 photometry 

The nucleus of M 31 is a photometrically and dynamically 
decoupled stellar system with respect to the surrounding 
bulge: it clearly stands out and shows no clear colour gra- 
dient from the UV to 1 /an except for a marginally resolved 
central UV excess interpreted as a compact stellar cluster 
(see e.g. Lauer et al. 1998 ). The detailed understanding 
of these observed structures first requires us to reconcile 
the miscellaneous published photometric and kinematical 
data sets. In this context, the WFPC2/HST images are key- 
stones in the study of the morphology of the nucleus of 
M 31. These data are now publicly available in a number 
of different bands. Photometric decomposition was done 
by Kormendy & Bender ( 1999[ hereafter KB99) and Lauer 
et al. (1998). We however decided to revisit the surface 



brightness decomposition to provide new quantitative es- 
timates regarding the respective contribution and colours 
of the bulge, the nucleus and the UV peak. An accurate 
centring of these structures is critical for the kinematical 
study conducted in the following Sections. 



2.1. The data 

We used WFPC2 images retrieved from the ST/ECF 
archives in the three following bandsQ: F300W, F555W 
and F814W (PI Westphal, 5236). The individual PC2 
exposures were combined via a drizzling technique 
(drizzle2) available under IRAF and implemented by H.- 
M. Adorf and R. Hook. This led to a pixel size of 0"022775, 
half a standard PC2 pixel. This had the advantage of 
correcting for the spatial distortions and partially solve 
the undersampling problem. The drizzling routine was ap- 
plied to all images using a new pixel fraction (set with the 
"pixfrac" parameter) of 0.6, except for the F300W expo- 
sure for which we kept a value of 1 as the individual PC2 
images were not dithered. The output Point Spread Func- 
tion (hereafter PSF) is roughly a quadratic sum of three 
contributions: the optical PSF, the input pixel size and 
the "pixfrac" parameter. The shorter wavelength of the 
F300W image thus compensates for the higher "pixfrac" 
value and gives an output PSF width comparable to the 
F555W drizzled image (about 1.6 PC2 pixels). All im- 
ages were carefully centred, normalised using the most 
recent PHOTFLAM values to convert the F300W, F555W 
and F814W into the VEGAMAG system. In the following we 
will simply write of the U , V and / images instead of the 
F300W, F555W and F814W respectively. We finally de- 
convolved these images with the corresponding PSFs (us- 
ing a Lucy-Richardson algorithm available in IRAF), and 
reconvolved them to a common resolution of about 0'.'06 
arcsecond, using a simple gaussian PSF. We checked that 
the residual differences in the final PSFs were small and 
do not affect the results presented below. 

2.2. Decomposition: the bulge, the nucleus and the UV 
peak 



As in Bacon et al. (1994), we used the Multi-Gaussian Ex- 
pansion formalism (Monnet et al. 



1992, Emsellem et al. 



1994) to build a two-dimensional fit of the surface bright- 
ness distribution in the centre of M 31. We first obtained 
an excellent fit of the bulge light in the / band using 3 
Gaussians. This model is only intended to fit the bulge 
surface brightness in the PC field. We also assumed a flat 
central brightness profile for the bulge by masking out the 
central 2" (where the nuclear light dominates) during the 
fitting process. The resulting parameters are presented in 
Table @. 

Our bulge representation differs from KB99's as the 
former has a non-zero and varying ellipticity, and a sur- 
face brightness profile which flattens towards the centre. 
The contribution of the bulge inside 2" is in fact impossi- 
ble to quantify given the present data and should there- 
fore be taken as an unknown free parameter of the model. 
We then normalised the / band bulge model to the U 



1 We also made use of the F568N exposures (PI Bohlin, 
5121): see Appendix A. 
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Table 1. Parameters for the MGE model of the bulge 
surface brightness in the PC field (/ band). 



# 


li 




1 


PA 




[L Q . P c- 2 ] 


[arcsec] 




[<ieg] 


1 


11012.11 


3.67 


1.000 


45.59 


2 


7709.34 


8.61 


0.667 


45.59 


3 


22970.38 


26.16 


0.964 


45.59 




-2 2 

x [arcsec] 



2 

y [arcsec] 



Fig. 1. Upper panels: major- (left) and minor-axis (right) 
V band surface brightness profiles with two models: our 
MGE fit (solid line), and KB99's spherical model (dotted 
line). The lower panels show the corresponding residuals 
after subtraction of the MGE bulge model. 



and V bands (see the comparison with KB99's model in 
Fig. |j]) and subtracted its contribution, providing bulge- 
subtracted images (see Fig. p|). As suming that the UV 
peak detected by King et al. ( 1995| ) has a negligible con- 
tribution in the / band, we normalised the nuclear / band 
image and subtracted it from the (bulge-subtracted) U 
and V images (Fig. [|). 

It is surprising to see how well this simple decomposi- 
tion procedure works: the only significant residuals in the 
U band indeed come from the UV peak (also seen in the 
V band) and from very blue (extreme horizontal branch?) 

1991 ). 



stars clustering around the nucleus (Brown et al. 
This shows that our bulge model is adequate (keeping in 
mind that the exact bulge light contribution remains un- 
known in the central arcsecond). It accounts for the ob- 
served ellipticity and position angle of the isophotes in 
the central 20 arcseconds. It also confirms that the core of 
M 31 can indeed be decomposed into three distinct com- 
ponents: the bulge, the nucleus, including PI, and the UV 




WFPC2/HST 
Bulge-subtracted 
nucleus of M 31 



Fig. 2. Bulge-subtracted isophotes of the nucleus of M 31 
from the WFPC2/HST images (top right: F814W /; top left: 
F555W V and bottom right: F300W U). The isophote step 
is 0.5 magnitude, and the faintest one is 15, 16.3 and 18.66 
for J, V and U respectively. The F300W faintest isophotes 
are significantly disturbed due to the lower S/N and the 
presence of blue stars. 

1.0 105 



60000 



60000 - 



40000 - 



20000 - 
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Fig. 3. Deconvolved bulge+nucleus subtracted surface 
brightness profiles of the so-called "UV peak" in the nu- 
cleus of M 31 (WFPC2/HST) in V (solid line) and U (dashed 
line). Flux units are in 2v0.pc -2 in the corresponding band 
(corrected for galactic extinction). 



peak. The obtained colours for the bulge and the nucleus 
are then (corrected for galactic extinction: 0.42, 0.24, and 
0.1436 mag in U, V and I respectively): (V — I)b = 1-26, 
(U-V) B = 1-98 and (V-I) N = 1.30, (U - V) N = 2.36. 
The UV peak has an average axis ratio of ~ 0.7 ± 0.05, a 
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PA of ~ 62° ± 8, a major-axis FWHM of 0"21, and has 
(U — V) = —0.35. We therefore find that the nucleus is 
redder than the bulge: this is in contradiction with the 
results of Lauer et al. ( 1998| ) who found that the nuclear 
region is bluer than the bulge in the U — V colour. We 
checked and confirm the robustness of this result by reex- 
amining the original images. We thus cannot explain this 
discrepancy even considering the slight differences in the 
reduction process. 



2.3. Centres and positioning 

Follo wing the photometric decomposition achieved in 
Sect. \i.2[ we have chosen our reference centre, position 
[0, 0] throughout the paper, to correspond to the central 
UV maximum in the F300W image (see Fig. ||). This 
choice was motivated by the fact that the UV peak is 
very bright and just resolved in the WFPC2 images, al- 
lowing a very accurate determination of its centre. It is 
also thought to correspond to the true position of the pre- 
sumed central massive black hole as discussed by KB99. 
KB99 defined their spatial zero radius x KB to be the ve- 
locity centre at the resolution of their SIS data. Compar- 
ing their surface brightness profile (their Fig. 8) with the 
F300W and F555W WFPC2 photometry we find that x KB 
is 0'.'031 away (towards PI) from the position of the UV 
peak (our [0,0]). KB99 quoted a value of 0'.'051 ± 0'.'014 
between the UV peak and x KB . This difference (of 0'.'02) is 
easily accounted by the fact that the contribution of the 
"UV peak" quickly decreases at longer wavelengths, the 
local maximum around P2 shifting away from PI, with an 
offset of about 0'.'023 already in the V band (roughly one 
dithered pixel). Throughout this paper, P2 is defined as 
the secondary surface brightness maximum in the / band, 
where we assume that the UV excess has a negligible con- 
tribution. P2 is indeed not coincident with the location 
of the UV peak: it shows up (in the / band) as an elon- 
gated structure, it centre being at about 0'.'076 from the 
UV peak (see Fig. ||). 



The spatial zero radius of Statler et al. (|1999| , Sta+99) 
is about — 0'.'025 from the UV peak (away from PI), ac- 
cording to their approximate surface brightness long-slit 
profile (spectral domain from ~ 4000 to 5450 A, see their 
Fig. 3). Their zero radius is therefore 0'.'025 + 0'.'031 = 
0'.'056 from the one defined in KB99. However, careful in- 
spection of Fig. 6 in Sta+99 shows that the kinematical 
data of KB99 has been uncorrectly shifted by ~ 0'.'13 (to- 
wards PI) in order to be consistent with their FOC kine- 
matics: this shift is then ~ 0'.'074 larger than it should be. 
The fact that Sta+99 managed to roughly fit the kinemat- 
ical data of KB99 comes partly from the way the extrapo- 
lation of their one-dimensional V and a FOC profiles was 
achieved, assuming almost no dependency perpendicular 
to the slit (see Sect. 5.1 of Sta+99), and partly from this 
uncorrect spatial registering. 



To summarize, the zero radius defined in Sta+99 and 
in KB99 are offset by -0"025 and 0'.'031, respectively, with 
respect to our spatial zero point taken as the location of 
the UV peak in the F300W band (Fig. |). We do not 
include here unknown (small) potential offsets perpendic- 
ular to their slits. 



3. OASIS spectrographic data 

3.1. Observations 

We observed the nucleus of M 31 in December 98 using the 
OASIS^ integral field spectrograph attached to the CFHT 
adaptive focus bonnette (PUEO). OASIS is the successor 
of the TIGER integral field spectrograph successfully used 
at CFHT between 1987 and 1996. The instrument design 



follows the TIGER concept described in Bacon et al. (1995), 
and includes a spectrographic mode, in which the spatial 
sampling is achieved via a microlens array, and an imaging 
mode. Various sampling sizes and spectral resolutions are 
available, from 0.04 arcsec to 0.16 arcsec per lens, and 
400 < R < 4000 respectively. OASIS has been designed 
and built at the Observatoire de Lyon and is operated at 
CFHT as a guest instrument. 

We used the MR3 configuration covering the Ca triplet 
(8500 A) region. This configuration was favoured against 
the classical Mg2 (5200 A) wavelength range because 
PUEO, like all other adaptive optics sys tems, performs bet- 
ter at longer wavelength (Rigaut et al. 1998 ). The selected 
spatial sampling of O'.'ll per (hexagonal) lens provides a 
field of view of 4x3 arcsec 2 . Details of the instrumental 
setup are given in Table ||. 

A total of nine exposures, each 30 mn long, were 
obtained during a seven day period (December 17-24). 
The atmospheric conditions were photometric. Seeing con- 
ditions were generally good, but changed rapidly (see 
Sect. 3.2.2). All exposures were centred on the nucleus, 
with only small offsets (typically of the order of 1-2 sam- 
pling size). Neon arc lamp exposures were obtained be- 
fore and after each object integration, and other required 
configurations exposures (bias, dome flatfield, micropupil) 
taken during daytime. Bright sky flatfields were also ob- 
served at dusk or dawn. The star HD 26162 (K2 III) was 
chosen as a kinematical template, and observed with the 
same instrumental setup. 

3.2. Data reduction 

The OASIS data were processed using the dedicated XOa- 
sis software (version 4.2) developed in LyonQ. The stan- 
dard OASIS reduction procedure includes CCD prepro- 



2 OASIS stands for Optically Adaptive System for Imaging 
Spectroscopy. It has been funded by the CNRS, the MENRT 
and the Region Rhone- Alpes 



3 Software and documentation are available at 


http: //www- 


]phs .univ-lyon 1 . fr / ~oasis /home /home_oasis_gb .ht 


mi 
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Table 2. OASIS instrumental setup 



PUEO 


Loop mode 




automatic 


Loop gain 




80 


Beam splitter 




I 


OASIS 



Spatial sampling 
Field of view 
Number of spectra 
Spectral sampling 
Instrumental broadening (a) 
Wavelength interval 



0.11 arcsec 
4x3 arcsec 2 
1123 

2.17 A pixel -1 
69 kms" 1 
8351-9147 A 



cessing (bias and dark subtraction), spectra extraction, 
wavelength calibration, spectro-spatial flatfielding, cosmic 
rays removal and flux calibration. Given the high surface 
brightness of the nucleus of M 31 and the small field of 
view, no sky subtraction was required. 

3.2.1. Spectra extraction and calibration 

The extraction of the spectra is the most critical phase of 
the reduction process. It uses a physical model of the in- 
strument with finely tuned optical parameters, such as the 
rotation angle of the grating with respect to the CCD col- 
umn, the tilt of the grism, the collimator and camera focal 
ratios, etc. The free parameters are fitted using a set of 
calibration exposures (micropupils - obtained without the 
grism -, arc and tungsten continuum frames) and saved in 
an extraction mask, later used to retrieve the precise (sub- 
pixel) location of each spectrum. Small shifts of the spec- 
tra location obtained at different airmass may occur due 
to mechanical flexures. These are corrected by including a 
global shift between the object exposure and the extrac- 
tion mask. This offset is measured via the cross-correlation 
of the arc frame associated with the science exposure, and 
the one associated with the reference continuum exposure 
(used to build the extraction mask). Shifts are typically 
0.1-0.2 pixel. Given that observations of M 31 were split 
over eight nights, we built two extraction masks from two 
sets of calibration exposures obtained at the start and end 
of the run respectively. These two masks were found to be 
almost identical, both giving an excellent fit (0.08 pixel 
RMS residuals between the model and the continuum ex- 
posure). Various extraction tests using each mask inde- 
pendently showed no significant differences, and we chose 
to use only the first mask. All spectra were extracted via 
an optimal extraction algorithm (Bacon et al. [2001 ) ap- 
plied to each exposure (object and associated arcs), finally 
providing 9 science datacubes with 1123 spectra each. 

These datacubes were then wavelength calibrated us- 
ing the associated (nearest) arc exposure. The goodness 
of the optical model guarantees that the raw extracted 
spectra do already have a good wavelength precalibration. 




-1 1 
x [arcsec] 



-10 1 
r [arcsec] 



Fig. 4. PSF fitting of OASIS exposure #6. Left panel: Con- 
tour plot of the reconstructed image (solid line) and the 
convolved HST image (dashed line). Right panel: photo- 
metric major-axis cut. 



First order correction polynomials were thus sufficient to 
finalise the calibration, with typical RMS residual errors 
of 0.04-0.05 A (~ 1.5 kms -1 ). After wavelength rebin- 
ning, the spectro-spatial flatfield correction was applied 
to each individual spectrum. This correction is obtained 
using a combination of a high S/N continuum datacube 
(tungsten lamp) and a twilight sky flat field datacube. The 
former is used to correct the spectral variations, while the 
latter allows to correct lens-to-lens flux variations. Cos- 
mic rays are then detected via a comparison beween each 
spectrum and its neighbours, and removed. The overall 
fraction of pixels affected by cosmic rays is small, typi- 
cally 0.2-0.3%. Finally, spectra are truncated to a common 
wavelength range (8361-8929 A) covering the Ca absorp- 
tion line triplet. No flux calibration was performed, since 
this is not required to measure the stellar kinematics. 

3.2.2. PSF determination, and merging of the datacubes 

PSF determination deserves special attention. Indeed, a 
precise knowledge of the spatial PSF is critical to be able 
to compare in detail datasets which span a rather large 
range in spatial resolution (from HST to ground-based 
long-slit observations). 

Reconstructed images are computed by direct summa- 
tion of spectra along the wavelength range followed by an 
interpolation on a square grid. As expected, the 9 recon- 
structed images present significantly different PSFs. We 
used the available HST/WFPC2 / band images of the M 31 
nucleus (see Sect. |J) to have an accurate estimate for each 
of them. Assuming a parametric shape for the PSF, we 
performed a simple least-squares fit between the convolved 
HST image and the equivalent reconstructed OASIS image. 
Free parameters for the fit also include the relative aligne- 
ment and rotation between the two images, as well as the 
corresponding flux normalisation factor. 

The reference HST/WFPC2 F814W (I) band filter covers 
a wavelength range which extends well beyond the spectral 
limits of the OASIS datacubes. However, colour gradients 
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Fig. 5. Fitted PSF of the nine OASIS exposures. FWHM 
range from 0779 (exposure #1) to 0'.'39 (exposure #6) 

Table 3. PSF parameters of the merged exposures 



ID 


Exp. 


FWHM 


o\ 


oi 


h/h 


0"3 


h/h 


M2 


6-7 


0.41 


0.15 


0.29 


0.98 


0.448 


0.023 


M8 


2-8 


0.50 


0.18 


0.35 


1.07 


0.407 


0.030 



in the M 31 nucleus are very small in the visible, so we 
decided to neglect the effects due to the differences in the 
response curve: we checked that this had a negligible effect 
on the final PSFs parameters. The PUE0 adaptive optics 
performance heavily depends on the "original" seeing con- 
ditions, the brightness of the guiding source and the wave- 
length range. OASIS observations of the M 31 nucleus are 
thus challenging for PUE0, given the low contrast and com- 
plex structure in the M 31 nucleus, and the relatively blue 
wavelength range of the observations (compared to the 
more classical near-infrared H or K bands). Typical PUE0 
PSFs at these wavelengths are not diffraction limited, but 
show a core of a few tenths of an arcsec surrounded by 
a large halo. A sum of 3 gaussian functions provides a 
reasonable approximation. An illustration of the quality 
of the fit is given in Fig. |4|. This procedure does not only 
provide PSF parameters but, as mentioned above, also the 
precise relative centring between the OASIS exposures and 
the HST/WFPC2 image. 

All nine computed PSFs are displayed in Fig. |H[ There 
are large differences in the resolutions of the nine expo- 
sures, which range from 0779 to 0/39 FWHM. Two merged 
exposures have thus been created using two different sets 
of exposures ranked according to their resolution: 

— the first one, hereafter called M2, with the highest spa- 
tial resolution is the combination of exposure #6 and 
#7 only. 



— the second one, hereafter called M8, with a higher S/N 
but a corresponding lower spatial resolution, is the 
combination of all exposures except exposure #1. 

The selected datacubes were merged using the accurate 
centring given by the PSF fitting procedure. Interpolation 
on a common square grid is performed before averaging 
the data. One lens located in the outer part of the nucleus 
was affected by a bad CCD column and was discarded be- 
fore interpolation. The PSF fitting was subsequently per- 
formed on the merged datacube and the results are given 
in Table [| The M2 and M8 exposures have cores FWHM 
of 0734 and 0741 respectively, both with extended wings 
giving a global FWHM of 0/41 and 0750 respectively. 



3.2.3. Computation of line-of-sight stellar velocity 
distributions 

A velocity template spectrum was obtained using the ref- 
erence star HD 26162 datacube, summed over an aperture 
of 1" to maximize the S/N. The spectrum was then contin- 
uum divided and rebinncd in ln(A). The same procedure 
was applied to the merged datacubes of M 31. We used 
the Fourier Correlation Q uotien t (FCQ) method, origi- 
nally developed by Bender ( 1990| ), to derive the kinemati- 
cal maps. The kinematics was parametrized using a simple 
Gaussian and complemented using third and fourt h orde r 
Gauss- Hermite moments (van der Marel & Franx 1993). 
All velocities were offset to heliocentric values. A systemic 
velocity of 308 km s _1 was deduced by comparing our data 
to the data of KB99 whose kinematics extends far enough 
to observe the slow rotation of the bulge (see Sect. ||). 



3.2.4. Bulge subtraction 

Since we are mostly interested in the kinematics of the 
nucleus of M 31, we would like to subtract the contri- 
bution of the bulge from the OASIS merged datacubes. 
This was performed using the photometric model derived 
in Sect. 2.2 (see KB99 for a similar approach). We ex- 
tracted a mean bulge spectrum from the outer part of the 
OASIS field, which we fitted making use of a library of 
stellar templates to get a high S/N reference spectrum. 
We then used a simple Jeans dynamical model (see Em- 
sellem et al. 1994) to derive the velocity and velocity dis- 



persion of the bulge (taking into account the instrumental 
setup and seeing) in the OASIS field of view. The resulting 
spectra were finally directly subtracted from the OASIS 
datacubes (M2 and M8) after the proper luminosity nor- 
malisation. We checked that slight variations in the contin- 
uum shape or absorption line depths did not produce any 
significant differences on the final bulge-subtracted dat- 
acube spectra. The results do however weakly depend on 
the details of the dispersion model used for the bulge: the 
same subtraction procedure was therefore also applied us- 
ing a constant dispersion of 150 krns" 1 for the bulge. The 
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and calibration files from the ST/ECF archives in this 
configuration^. After the data reduction provided via the 
ST/ECF pipeline using the best reference files, we cor- 
rected the individual frames for the residual fringing us- 
ing the dedicated routines available under IRAF (Goud- 
frooij & Christensen 1998) and the available contempo- 
raneous flat fields. We combined the individual exposures 
after careful recentring, and extracted individual spectra, 
which were rebinned in In (A) for the kinematical measure- 
ments. The spatial centre of the long-slit luminosity profile 
was then accurately determined taking the same reference 
zero point as in the WFPC2/HST data (see Sect. |]|). Spec- 
tra were binned along the slit to increase the S/N outside 
0'.'85 from the centre. We also retrieved STIS exposures 
from the K0III star HR 7615 (same configuration) from 
the archive to be used as a stellar template, and reduced 
them in the same way. The line-of-sight velocity distri- 
bution at each spatial location was finally derived using 
the Fourier Correlation Quotient method as in Sect. |3.2.3j 
We estimated a maximum velocity error of ~ 8 kms -1 
due to the slit effect using the present STIS characteris- 
tics (see Bacon et al. 1995| ). In the present paper, we will 
mainly focus on the first two moments, namely the ve- 
locity and velocity dispersion, only using the higher order 
Gauss-Hermite moments to derive an estimate of the first 
two true velocity moments. We finally performed the sub- 



traction of the contribution of the bulge as in Sect. 3.2.4 



Fig. 6. The central region of the deconvolved WFPC2 im- 
age in the I band showing the location of PI, P2 and 
the UV peak. Isophotes are drawn from 11.79 to 12.49 
with a step of 0.05 magarcsec -2 . The UV peak, located 
at [0,0], is marked by a cross. The solid line represents 
the kinematical major-axis, with which P2 elongation is 
well aligned. The slit axis for the SIS (KB99), STIS (this 
paper) and FOC (Sta+99) data are represented by dotted- 
dashed, dashed and dotted lines, respectively. The filled 
and empty circles represent the location of the spatial zero 
point of the SIS and FOC long-slit data, respectively. The 
upper panel shows a zoom in the central 0'.'2. Note the 
distance between the maximum of PI (marked with the 
text "PI") and the kinematical major-axis. 



difference comes mainly from the higher assigned ce ntral 
dispersion to the bulge in the Jeans model. In Sect. 5.2 , 
we will only discuss results from the bulge-subtracted M8 
OASIS datacube, which has a significantly better S/N. 



4. STIS/HST data 

The nucleus of M 31 was observed with STIS at PA= 39° 
in July 1999, using the G750M grating and the 52 x O'.'l 
slit (Proposal 8018, PI Green). The spatial pixel was 
0'.'05071, and the velocity resolution about 38 kms -1 . 
We have retrieved the 8 corresponding science exposures 



5. Results 

5.1. OASIS stellar velocity maps 

In Fig. [7|, we present the stellar velocity and velocity dis- 
persion maps^j for the two final datasets M2 and M8. Both 
give similar results, as shown in Fig. ||, although the dif- 
ference in spatial resolution can be seen in the central 
velocity gradient and velocity dispersion peak. 

We compare these results with the KB99 kinemati- 
cal profiles obtained with the SIS/CFHT spectrograph at 
a spatial resolution of 0'.'65 FWHM. The SIS slit posi- 
tion has been spatially shifted by +0'/031 to account for 
the small offset between the location of the UV peak and 
the velocity centre as defined by KB99 (see Sect. 2.3). 
A simulated OASIS slit was computed using the surface 
brightness weighted average of the kinematical quantities 
(V and yV 2 + a 2 ) over an equivalent slit opening of 0'.'35 
at PA= 52.5° (Fig. |9[). The comparison was done using 
the M8 data set. The agreement is excellent, with 8 and 
9 kms -1 RMS residual in velocity and velocity dispersion, 
respectively. Note that the difference due to spatial reso- 
lution is largely smoothed by the averaging over the 0'.'35 
slit width. 



4 We did not include here the G430L exposures as we were 
solely interested in the kinematics. 

□ASIS data are available in ascii/fits form at 



tittp : //www-obs .univ-lyonl . f r/~tigerteam/bacon01 .html 
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Fig. 7. OASIS stellar kinematical maps: stellar velocities (left panels) and dispersions (right panels). The top and 
bottom panels correspond to the high spatial resolution datacube (M2) and the high S/N datacube (M8), respectively. 
The dotted white lines display the zero isovelocity contour and the 200 kms -1 isovelocity dispersion. The step is 
25 kms -1 in all panels. The centre of the UV peak, our [0, 0], is marked by a filled circle, the cross marking the centre 
of PI. The dashed lines indicate the kinematical major-axis. North is at a PA of 192°, East at its right, and the scale 
is in arcsecond. 



The precise axis and centre of symmetry of the ve- 
locity field was estimated using a fit with a simple an- 
alytical function. In that model, the velocity field was 
parametrized as a sum of first order Gauss-Hermite-like 
functions x' ■ J2 i Gauss^x', y') where x',y' are cartesian 
coordinates on the sky. Free variables are the coordinates 
of the individual centres (xq, j/q), the principal axis angle 
(#o)i the Gaussian parameters (o-j, Ij, and a zero veloc- 
ity. Experiments show that the kinematic centre and tilt 
are insensitive to the details of the fitting functions pro- 
vided it is antisymmetric. The two merged datacubes gives 
consistent results. A use of only two components provides 
a reasonable fit with 16 kms -1 RMS residuals. The cen- 
tre of symmetry is found at [0,0] within 20 mas and the 
kinematic axis with 9q = 56.4 ± 0.2°. The velocity at [0,0] 
is not zero, but ~ —9 kms -1 . 

The kinematic axis is significantly different from the 
P1-P2 axis (PA 42°) as shown in Fig. g, and PI is off- 
set 0'.'12 from the kinematic axis. The fact that PI is not 
aligned on the kinematic major-axis must be taken into 
consideration for the interpretation of the high spatial res- 



olution HST kinematical data which have been taken close 
to the P1-P2 axis. 



The velocity dispersion map of the (best resolved) M2 
datacube peaks at 270 kms -1 at 0'.'3 ± O'.'l on the anti- 
Pi side. The peak is extended, and surrounded by a halo 
which itself extends above the nearly constant bulge ve- 
locity dispersion of ~ 150 kms -1 . The same structure is 
found in the M8 datacube, but with a lower contrast. The 
dispersion is clearly asymmetric with respect to the UV 
peak, with a difference of 35 kms -1 at ±l'/2 along the 
kinematic axis. Note that the offset of the velocity dis- 
persion peak is also present in the KB99 data (see their 
Fig. 4). The small difference is simply due to the slightly 
lower resolution of the SIS data (including the smoothing 
onto the 0'.'35 arcsec slit width). Such an offset was also 



present in the TIGER data (Bacon et al. 1994 ), but with a 
larger value (0'.'7), due to the lower spatial resolution (0'.'9 
FWHM). 
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Fig. 8. OASIS velocity profiles along the kinematic axis 
(PA= 56.4°, averaged over a width of 0'.'2). Top panel: 
stellar velocity. Bottom panel: stellar velocity dispersion. 
The crosses and circles correspond to the high S/N (M8), 
and the high spatial resolution (M2) datacubes respec- 
tively. 



5.2. Bulge- subtracted velocity maps 

We have used two different models for the (unknown) 
velocity dispersion of the bulge: a constant value of 
150 kms -1 , or the dispersion predicted by a simple Jeans 
model using the combined gravitational poten tial of the 
nucleus and the bulge (as in Bacon et al. 1994 ). The two 
resulting kincmatical profiles only show minor differences, 
the bulge-subtracted velocity and dispersion for the "con- 
stant dispersion" model being slightly larger: for clarity, 
we will only deal with the latter (note that these pro- 
files are consistent with the bulge-subtracted kinematics 
of KB99). 

After the subtraction, the asymmetry in the veloc- 
ity field is now clear, with the anti-Pi side having a 
slightly higher velocity amplitude with a local minimum 
at -250 kms -1 (at - 0"75 from the centre; fig. [l(]). The 
velocity profile on the PI side is nearly flat with a value of 
~ 180 kms -1 . The dispersion now peaks at 329 kms -1 at 
— O'.'l from the centre. We also confirms that the nucleus 
is cold (KB99), with a value of 108 ± 10 kms -1 at l'/2 
along the kinematical major-axis on the PI side. 



> 



100 


-100 
250 



g 200 
b 150 h 



100 



_ 1 1 1 1 1 1 1 


+ 
m 

+ 


i i I i i i _ 

+ * _ 


!» + + 

" ++ +„ + 


— 1— r 


PI - 


; I I I | I I I 
: t 

_ + 

: + ++++ ++* 
:i r 


f + 

i 


— 1 1 1 1 1 1 — \ 

+i V+ s + * 7 
+ » +* + : 

_ 



x [arcsec] 

Fig. 9. Comparison between the SIS kinematics (circles; 
Kormendy & Bender |1999| ; PA= 52.5°, slit width of 0"35) 
and the OASIS equivalent slit (crosses). Top panel: mean 
velocity. Bottom panel: velocity dispersion. 

5.3. STIS and F0C kinematical profiles 

The STIS velocity and dispersion profiles^ are presented in 
Fig. [ll]. There is a clear asymmetry in the velocity curve, 
with local turnover values of 197±5 and — 292±20 kms -1 
at radii of +0'/42 and — 0'.'235 respectively. The maximum 
velocity on the PI side is 201 ± 5 kms" 1 at 0"63. The 
dispersion is maximum on the anti-Pi side at a radius of 
— 0'.'235 with a value of 321 ± 33 kms^ 1 (and a value of 
313±36 kms -1 at — 0//18). This is to be compared with the 
offset of 0'.'15-0'.'20 founcQ by KB99 and - 0'.'2 for the slit 
profile reconstructed from the OASIS data. The STIS ve- 
locity profile crosses the V = line at +0//09 (PI side). At 
[0,0] we measure V = — 68 ± 6 kms -1 . After subtraction 
from the bulge contribution, the central velocity gradient 
is slightly steeper, with turnover values of —355 kms -1 
and 221 kms" 1 at -0'.'235 and 0"47 respectively. 

Comparing the original F0C (Sta+99) and the STIS 
data, we find significant discrepancies in both the veloc- 
ity and dispersion profiles at the very centre which seem 
difficult to attribute to differences in instrumental char- 
acteristics. The most surprising difference is the location 
of the dispersion peak in the original F0C data which is 
nearly centred with an offset of only 0'.'06 from the UV 
peak (away from PI). By including the observed spatial 
shift of 0'.'025 mentioned in Sect. |2.3] , as well as a veloc- 
ity shift of 30 kms -1 , the comparison looks reasonable 



(Fig. |l2j): this point is discussed below (Sect. 5.4). 
STIS data are available in ascii/fits form 



at 



[iittp : //www-obs .univ-lyonl . f r/^tigerteam/bacon01 .html 

7 T his includes the spatial offset of / /031 discussed in 
Sect. El 
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Fig. 10. Bulge-subtracted OASIS kinematics along the 
kinematic major-axis of the nucleus compared with the 
original OASIS M8 profiles (crosses). The solid line shows 
the resulting kinematics when using a constant dispersion 
of 150 km s _1 for the bulge, and the dotted line when using 
the dispersion predicted by a simple Jeans model. 



5.4- OASIS compared with HST kinematics 

Although the STIS and FOC data have similar slit widths 
(O'.'l for STIS and 0'.'063 for FOC) and PAs (39° for STIS 
and 42° for FOC), the FOC data (Sta+99) have been ob- 
tained at much bluer wavelengths than STIS (0.45 versus 
0.85 fj,m) and have a better spatial resolution (Fig. 12). A 
detailed comparison between the STIS and FOC data sets 
thus requires to take into account these specific charac- 
teristics. But this implies the convolution of the (yet un- 
known) 2D kinematics of M 31 at FOC spatial resolution. 
We also wish here to compare HST and OASIS kinematical 
data: the two-dimensional coverage of the OASIS spectra is 
in this context a powerful additional constraint. The com- 
parisons between the OASIS and STIS datasets are shown 
in Fig. [L3|. 

To attempt such a comparison, we have built simple 
2D parametric kinematical models. These models are ob- 
viously ill-constrained, given that we have knowledge on 
only a tiny fraction (2.5%) of the area covered by OASIS at 
FOC resolution, and do not stand on any physical ground. 
This ad hoc model is solely designed to check if the three 
kinematical data sets are consistent with each other, i.e. if 
we can find a reasonable model which simultaneously fits 
all the observed kinematics. All convolutions are achieved 
taking into account the relevant pixel integration. The sur- 
face brightness distributions were derived from the decon- 
volved WFPC2 images: we kept the F814W filter image for 



200 




x [arcsec] 

Fig. 11. STIS kinematical profiles: velocity (top panels) 
and dispersion (bottom) panels at a PA of 39°. The in- 
serted panels present the STIS kinematical profiles within 
a radius of 10". 



both the STIS and OASIS data, and used a weighted sum 
of the F300W and F555W images to roughly reproduce 
the surface brightness profile of the FOC data presented by 
Sta+99. In the cases of the OASIS and STIS data, we used 
an estimate of the first two true velocity moments (the 
best gaussian velocity and velocity dispersion depending 
on the shape of the line-of-sight velocity distribution). For 
the FOC data, we had to use the original values as pub- 
lished by Sta+99. Since Sta+99 mention that the kine- 
matics derived from the half blue and red parts of their 
spectra are indistinguishable, wc assumed in the following 
that the UV peak does not contribute to the kinematical 
profile. In other words, we assumed that the stars form- 
ing the UV peak do not show (metallic) features strong 
enough to significantly contribute to the observed kine- 
matics: this hypothesis seems to be supported by a prelim- 
inary analysis of the STIS/G430L data (Emsellem et al., 
in preparation). 

We start with a model based on velocity functions 
of a simple form V(r) = K ■ x ■ (1 + r 2 /r 2 )~ n with 
r 2 = x 1 + y 2 /q 2 (x and y being the cartesian coordinates 
aligned with the kinematical axis), and constant or gaus- 
sian functions for the velocity dispersion. The PA of the 
kinematical major-axis was set to 56.4° as measured from 
the OASIS data. The observed asymmetries in the kinemat- 
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Fig. 14. Comparison of OASIS (left), STIS (middle) and FOC (right) kinematics within the central arcsec. Velocity and 
dispersion profiles are shown in the bottom and top panels respectively. Original FOC data (right panels) have been 
shifted by 0'.'025 and 30 kms" 1 (see text). The solid lines represent the model convolved and sampled according to 
the corresponding instrumental setup. 



ical profiles were obtained by allowing both spatial and ve- 
locity offsets in the parametrized functions. We fixed the 
dispersion of the bulge and the nucleus to 150 kms -1 and 
108 kms -1 respectively, adding (quadratically) a gaussian 
of l'/2 FWHM with a maximum of 140 kms -1 located at 
0'.'8 on the PI side to reproduce the low- frequency spa- 
tial asymmetry seen in the OASIS dispersion map (Fig. fjj). 
The velocity and the non-centred second order moment 
(V 2 + a 2 ) are then convolved, after proper luminosity 
weighting. The input parameters were tuned until a satis- 
factory result was obtained. 

The best overall fit was obtained by adding high veloc- 
ity components in the central 0'.'3 aligned with the kine- 
matic major-axis. The dispersion peak then simply results 
from the effect of velocity broadening as shown in Fig. |l5| . 
We cannot distinguish between a velocity broadening ef- 
fect from true velocity dispersion of components that are 
spatially unresolved. This model is anyway the simplest 
we found which provides a reasonable fit to the STIS and 
OASIS data. It includes maximum amplitude velocities of 

~ +235 and 265 kms" 1 (outside the central 0'.'3). We 

also confirm the spatial and velocity offsets in the origi- 
nal FOC data with respect to the UV peak and the zero 
velocity as measured by KB99, respectively: by shifting 
the FOC kine mat ics by 0'.'025, consistently with the anal- 
ysis of Sect. |2.3| , and adding 30 kms -1 to their systemic 



velocity, we can indeed reconcile the model with the FOC 
data (Fig. |l4|). Note that this velocity shift is consistent 
with the maximum possible systematic error of 50 kms -1 
quoted by Sta+99. The high velocity component of the 
model is an attempt to reproduce the abrupt jump seen 
in the FOC rotation curve at a radius ~ — O'.'l, near the 
location of the dispersion peak. This feature should how- 
ever not be overinterpreted, although we believe that the 
central dispersion peak can indeed be explained via veloc- 
ity broadening with the contribution of fast moving stars 
within 0'.'3 of the UV peak on the anti-Pi side. We do not 
find any significant discrepancy in the dispersion profiles 
contrarily with what was advocated by Sta+99: as men- 
tioned earlier (Sect. ^), they uncorrectly shifted KB99's 
data to manage a consistent fit (their Fig. 6). 

We cannot make any definite statement about the de- 
tailed kinematics in the central 0'.'3, as there is too much 
freedom in the model parameters. This is mainly due to 
the fact that FOC and STIS provide only one-dimensional 
profiles, unlike OASIS data which are two-dimensional, but 
have a too low spatial resolution to constrain the dynamics 
at that scale. We can however still make a few relatively 
safe statements, which will be discussed further in Sect. ^| 

— The observed kinematics is definitely not consistent 
with a dispersion peak at the location of the UV peak 
(our [0,0]). 
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Fig. 12. Comparison between STIS (solid line) and FDC 
kinematics (filled triangles). A shift of 0'.'025 towards the 
anti-Pi side was applied, following the analysis made in 
Sect. 2^. A positive shift of 30 kms" 1 was applied to the 
F0C velocity profile (top panel, see text) 
. The STIS data and the F0C data were taken at PAs of 39° 
and 42° respectively. 

— The dispersion peak and high velocities seem to be 
spatially associated with P2 (as seen in the I band) 
which is offset from the UV peak on the anti-Pi side 
as shown in Fig. [| 

- The OASIS, STIS and F0C kinematical profiles are con- 
sistent with each other, although a two-dimensional 
coverage of the central 0'.'3 is required to properly ad- 
dress the kinematics there. 

6. Numerical Simulations 

We now interpret the observations with the help of N- 
body modeling. Among the various hypothesis that have 
been advanced to explain the M 31 double nucleus (with 
luminosity peak PI shifted by ~ 1.8 pc from the kine- 
matical centre, almost coinciding with the UV peak), the 
most natural would be an m = 1 wave in a rather cold 
and thin disk orbiting the SBH, located at the centre of 
the UV peak. The disk has to be cold, to be unstable to 
non-axisymmetric waves, and this implies a small thick- 
ness, in the almost spherical potential provided by the 
central SBH. However it is still unclear whether m = 1 
perturbations, accompanied by a displacement of the grav- 
ity centre, will be unstable, and develop spontaneously in 
the physical condi tions corresponding t o the M 31 nucleus 
(Heemskerk et al. 1992 , Lovelace et al 1999| ). 
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Fig. 13. Comparison between the kinematics from STIS 
(crosses) and OASIS (filled circles). The OASIS kinematics 
have been averaged over a 0'.'2 wide slit (PA= 39°). 




Fig. 15. Same as in Fig. |14| but for the OASIS and STIS 

data in the central 2". 



An alternative solution is that the disk is initially per- 
turbed by an external force (either a globular cluster or 
Giant Molecular Cloud passing by), and the response is 
long-lived with a time-scale comparable to that of such 
external perturbations. We have simulated this possibil- 
ity, and indeed found modes of m = 1 oscillations that 
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maintained during 70 Myr, with an almost constant pat- 
tern speed. This peculiar feature is due to the very low 
precession rate in an almost Keplerian disk, near a SBH. 
The asymmetry of the density, and the radial variation 
of eccentricity of the orbits, generate local variation of 
the effective precession rate, and most of the stars are 
dragged into a mode of slow, positive pattern speed. In 
this paper, we propose that this mechanism is at the ori- 
gin of the M 31 m = 1 perturbation (or "double-nucleus" 
morphology), and illustrate this with N-body simulations 
with asymmetric initial conditions. We consider in a fu- 
ture work the possibility that the m = 1 instability de- 
velop spontaneously from axisymmetric initial conditions 
(Combes & Emsellem |200l| ). 



6.1. N-body Methods and Diagnostics 

We performed essentially 2D simulations, since our inter- 
pretation is that the nuclear disk of M 31 is cold and 
thin (axis ratio of the order 0.1). The apparent axis ra- 
tio must then be mostly due to an inclination of i ~ 55° 
(see Sect. |6^ ), less edge-on than the large-scale M 31 disk 
(i = 77°). However we also checked the results with 3D 
experiments. The gravity is solved via fast Fourier Trans- 
forms, on a useful 2D grid of 256 x 256 (512 x 512 to sup- 
press Fourier images). In 3D, the N -body code is a FFT 



scheme, which uses the James (y_977j) method to avoid the 
influence of Fourier images. The grid is then 128 x 128 x 64. 
The size of the simulated box in the disk plane is 20 pc, 
corresponding to a cell size of 0.078 pc (and twice that in 
3D). This size is also the softening length of the Newtonian 
gravity. The M 31 nuclear stellar disk is represented in 2D 
by 99000 particles and 152384 in 3D. Two rigid potentials 
arc added, representing the bulge and the supermassive 
black hole. The time step is 10~ 4 Myr. 

Particle plots are made in face-on and M 31-sky projec- 
tions, together with the velocity field, the density, velocity 
and dispersion profiles along the major axis, to compare 
with the observed quantities (with and without bulge ad- 
dition). The Fourier analysis of the density and potential 
are made regularly as a function of time and radius. If 
the potential is decomposed as 3>(r,0) = $o(r) + < I ) m( r ) 
cos (m 8 - 4> m ), we define the intensities of the various 
components by their maximal contribution to the tangen- 
tial force, normalised by the radial force F r = —d&o/dr, 
through 

S m = m$ m /r\F r \ 

In some simulations the SBH is allowed to move, with 
respect to the fixed background bulge. However, this in- 
troduces additional oscillations, of small amplitude, which 
are not well computed with the low spatial resolution in 
3D. Therefore in most cases, the SBH was fixed at the cen- 
tre of the bulge. These oscillations are of different (higher) 
frequencies that the phenomenon we are studying here, 
and are only superposed to it. They may be related to 



the nuclear oscillations studied by Taga & lye ( 1995 ) and 
Miller & Smith fll99"^ ). 

Over all run periods of typically 10-100 Myr, the 
Fourier analysis of the disk potential was done as a func- 
tion of radius, and stored every 0.03 Myr. The result 
F(r, t) is then Fourier transformed in the time dimension, 
in order to get the power as a function of pattern speed, 
F(r, ft). With respect to the frequency of the m perturba- 
tion, the pattern speed is O = uj/m. 

6.2. Galaxy Model and Initial Conditions 

The massive black hole potential was softened to a few 
cell sizes to avoid prohibitively small values of the time 
step in the centre. It is represented by a Plummer shape 
potential, 

GM B h 



V r2 + ■ 



bh 



with ri,h = 0.07 pc. The mass of the SBH Mbh was varied 
between 3.5 and 10 10 7 M Q to test the role of the self- 
gravity of the disk. The fixed analytical potential of the 
bulge is composed of 4 Plummer functions, determined 
from the MGE method (see Emsellem & Combes 1997). 
Their corresponding masses and radii are displayed in Ta- 
ble £| The total mass of the bulge inside 9 pc is 0.87 10 7 
M . The nuclear stellar disk is initially a Kuzmin-Toomre 
disk of surface density 

£(r) = £ (l + r 2 /d 2 r 3/2 

truncated at 9 pc, with a mass of Mj, = 1.7 10 7 M Q , and 
characteristic radius of d = 3.5 pc. It is initially quite cold, 
with a Toomre Q parameter of 0.3. The mass fraction of 
the disk inside 9 pc was varied between 20 and 40%. 

Since the potential of the black hole is significantly 
smoothed below a radius equal to 2.5 times the soften- 
ing length, here 2.5 x 0.07 pc= 0.175 pc, we have avoided 
putting particles in this region, so as not to introduce arti- 
ficial dynamics. The initial surface density was then equal 
to the difference between two Toomre disks of the same 
central surface density, but different characteristic lengths. 
This was done to have an analytic density-potential pair 
with a smooth boundary for the central hole. Up to 12% of 
the central disk mass was thus removed and redistributed 
on the remaining outer parts. The exact value of the scale- 
length of the removed disk, between 0.5 and 1.2 pc, was of 
no importance on the final results. It is somewhat larger 
than the required 2.5 softening lengths, since the initial 
orbits near the black hole are eccentric, with a small peri- 
centre. The potential inside this radius is in any case dom- 
inated by the keplerian law of the black hole, so that the 
suppression of the central stellar disk does not affect the 
rotation curve. The rotation curve obtained is plotted in 

Fig. m 

The initial m = 1 perturbation was produced in several 
ways: either the SBH position was shifted by an arbitrary 
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Table 4. Parameters of the mass model 



6.3. Results 



Component 


Size 


Mass* 




pc 


10 7 M Q 


SBH 


0.07 


7 


Disk 


3.5 


1.7 


Bulge- 1 


700 


0.01 


Bulge-2 


160 


0.07 


Bulge- 3 


50 


0.2 


Bulge-4 


15 


0.6 



The best fit model for M 31 was obtained through the 
initial distribution of eccentricity 



mass inside 10 pc radius 
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Fig. 16. The observed orbital velocity V (full line, in 
km/s), log of surface density in arbitrary units (S, dash), 
and radial velocity dispersion (a v , dot-dash), at the epoch 
of 54 Myr, for the galaxy model with M B h = 7 10 7 M . 



value from the centre of the potential, or its velocity was 
perturbed and given a high value. In both long- 
lived m — 1 pattern was obtained. However, the initial 
conditions were too far from equilibrium, and the violent 
relaxation led too many particles to escape. A softer way 
to perturb the disk is to launch particles in eccentric orbits 
from the start, keeping the SBH fixed at one of the foci of 
the elliptical orbits. Once a profile of eccentricity is chosen, 
each particle is moved along the corresponding ellipse at a 
random position with probability inversely proportional to 
its linear velocity at this point (probability density along 
the ellipse proportional to 1 — ecos(£), £ being the true 
anomaly). About 20 models were run, to investigate the 
nature of the modes, and the influence of initial conditions, 
the mass fraction of the disk, etc. The best fit for M 31 is 
described here, and compared with observations. 



e(a) = 0.6^/(1 -(a/4 pc) 2 ) 

and e(a) = for a > Ape (a being the semi- major axis). 
The major axis of the orbits are initially aligned, but they 
quickly re-arrange, so that the pericentre of orbits in the 
inner and outer parts are in phase opposition. The face-on 
view of the nuclear disk is displayed in Fig. |l8|. There is 
a clear density accumulation on one side of the BH. The 
strength of the m — 1 and m = 2 Fourier components 
of the potential are plotted in Fig. |l^, together with the 
position of the stellar disk gravity centre, as a function 
of time. At the beginning, until 15 Myr, there is a set- 
tling phase, where the m — 1 mode has a trailing, than 
clearly leading arm. The stellar gravity centre precesses in 
the positive direction, and also the pattern speed in the 
power spectrum is positive (see Fig. |T^). Then, the motion 
of the gravity centre becomes more regular, slows down 
slightly, with a period of 2 Myr, corresponding to a pat- 
tern speed of the m = 1 perturbation of fl p ~ 3 km/s/pc 
(the pattern speed is equal to the orbital frequency of the 
gravity centre). Note that the m = 2 perturbation is just 
the harmonic of the m = 1, and has the same pattern 
speed. 

These values of pattern speed are about 10 times 
smaller than the value derived by Sambhus & Sridhar 
(2000) using the Tremaine- Weinberg method and M 31 
FOC data. Note however that the latter method is rather 
uncertain when applied to one-dimensional data. 

The distribution of eccentricity is quite similar all over 
the run, and typically represented by Fig. [2^. It is always 
decreasing from a = 1 pc to 6 pc, with a large scatter. The 
eccentricity is maintained at large values, and the apoc- 
cntres are aligned to give the density accumulation that 
will give rise to PI in M 31 (see the position of PI in the 
middle and bottom of Fig. g^). From 5 to 10 pc, the apoc- 
entres change side, and there is a secondary density ac- 
cumulation, in phase opposition with SI, of much smaller 
amplitude. The accumulation is also at the apocentre of 
these outer particles. The density accumulation is located 
very close to the maximum eccentricity gradient. As noted 
by Statler (1999), this is a configuration that allows the 
right impulses on eccentric orbit to equalize the precessing 
rates. If an orbit grazes from the inside a density accu- 
mulation, the experienced outward pull can accelerate in 
the positive sense its precession, if it is at apocentre. If it 
is at pericentre, its precession will be slowed down. The 
nearly axisymmetric (epicyclic approximation) precessing 
rate is only slightly negative, and this small impulse is suf- 
ficient to align the orbits, and make them precess in the 
positive sense. The orbits should be at their apocentre in- 
side the density peak, and at their pericentre inside, for 
the precessing corrections to go towards equalization. This 
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Fig. 17. Pattern speed as a function of radius, in units of 
km/s/pc, for the m = 1 mode (top) and the m = 2 mode 
(middle) between the epochs 28.8 and 43.2 Myr (Mbh = 
7 10 7 M ). The pattern speed slows down slightly with 
time, which thickens the lines (from 3 to 2.5 km/s/pc in 
60 Myr). 




-5 







Fig. 18. Face-on surface density of the nuclear stellar disk, 
at epoch 28.8 Myr (M BH = 7 10 7 M Q ). The scale is in 
pc. The hole in the centre is introduced on purpose in 
the initial conditions, to avoid artificial effects due to the 
black hole softening and limited spatial resolution (see 
Sect, f 



explains the orbit orientation observed, and the fact that 
the two density accumulations are in phase opposition. 

The morphology and kinematics of the simulated nu- 
clear disk at epoch T — 10 Myr is compared to the M 31 
observations in Figs ||FJ and ^l|. With an inclination on the 
plane of the sky of i — 55° ± 5, and a position angle of the 
density maximum such that PI is not on the line of nodes, 
the comparison is reasonably satisfying. The constraint on 
the inclination angle is obviously model dependent. How- 
ever, the disk needs to be thin for non-axisymmetric wave 
to develop. For inclinations i > 60°, it will thus be hard 
to reconcile our proposed model with the observations. 
There is of course a deep hole in the centre of the simu- 
lated image, since we have avoided placing particles here 



(see Sect. 6.2). The limited spatial resolution imposes a 
softening of the BH to 0.1 pc scale, and prevents follow- 
ing the stellar dynamics there. For Figs |2TJ and [HI we 
have artificially filled in this hole by adding a central ax- 
isymmetric disk of particles: they do not participate in 
the m = 1 mode but their velocities are of course chosen 
according to the background potential. 

Although the presented model has many limitations, it 
reproduces the most important observed features, specif- 
ically the asymmetries in the surface brightness between 
PI and its phase opposite (Fig. |2^) as well as in the ve- 
locity profile (Fig. 21). A thorough comparison, e.g. in- 



important comment concerning stars in the central par- 
sec. We expect the stars in this region to participate to 
the m = 1 mode, as observed in self-consistent simula- 
tions where we initially did not remove the central parti- 
cles. In that case the sign of the eccentricity was reversed 
in the central parsec, creating an offset light peak on the 
anti-Pi side. This configuration was already suggested by 
Statler ( 1999| ) who studied the sequence of closed orbits 
in a wide precessing disk. It is consistent with the ob- 
served offsets of P2 (/ Band) surface brightness, and of 
the dispersion peak on the anti-Pi side (see Sect. 5.4 and 
Fig. |J). The UV peak would thus still mark the location 
of the SBH, but it would not contribute to the measured 
kinematics. 



6.4- Discussion 

The m = 1 modes of a self-gravitating disk, dominated 
by a central point mass, have not been fully investigated. 



eluding the central velocity dispersion profile, should wait 
for simulations including a more realistic treatment of the 
particles near the SBH. We can however already add an 



Jacobs & Sellwood (1999) have reported the presence of 
a slowly decaying m = 1 mode in annular disks around 
a slightly softened point mass, but only for disk masses 
less than 10% of the central mass concentration. It is ex- 
pected that there is a threshold mass of the central BH, 
above which eccentric modes could be sustained by grav- 
ity. Indeed, in the extreme cases, where the disk mass is 
negligible, the orbits are Keplerian, and their precessing 
rate is exactly zero (fi = k). If the apsides are aligned at 
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Fig. 19. top: Coordinates (X: full line) of the centre of 
gravity of the stellar disk, as a function of time (model 
with Mbh = 710 7 M Q ). Note that the period is slightly 
increasing with time, i.e. the pattern speed slows down, 
bottom: Intensity of the to — 1 (SI, solid line) and to = 2 
(S2, light line) Fourier components of the potential (more 
exactly the corresp. components of the tangential force 
normalised by the radial force). 



a given time, they will stay so in a Sl p = mode. The 
self-gravity of the disk makes k > SI, and the orbits dif- 
ferentially precess at a rate SI — k < . However, if the 
disk self-gravity is not large, a small density perturbation 
could be sufficient to counteract the small differential pre- 
cession. Goldreich & Tremaine ( |1 9 79 ) showed that in the 
case of Uranian rings, the self-gravity could provide the 
slight impulse to equalize the precessing rates, and align 
the apsides. Levine & Sparke ( 1995 ) proposed that lopsid- 
edness could survive if the disk is orbiting in an extended 
dark halo, provided it remains in the region of constant 
density of the halo (or consta nt SI, but this does not ap- 
ply to M 31). Lovelace et al. ( 1999| ) found through linear 
analysis slowly growing modes, in the outer parts of a disk 
orbiting a central point mass. The patter n spee d could be 
either positive or negative. Taga & lye (1998) found by 
N-body simulations that a massive central body can un- 
dertake long-lasting oscillations, but only when its mass 
is lower than 10% of the disk mass (which again does not 
apply to M 31). 

Here we propose that a natural m = 1 mode can ex- 
plain the M 31 eccentric nuclear disk. We have found that 
for a disk mass accounting for ~ 20 - 40% of the total cen- 





Fig. 20. Comparison between the bulge-subtracted 
WFPC2/HST / band surface brigtness of the nucleus of 
M 31 (thin solid line) and an N-body simulation of an 
to = 1 mode where the disk was 40% of the central mass 
(with and without the addition of particles in the central 
hole as mentioned in the text, represented by dotted and 
thick solid lines respectively). The inclination of the disk 
was chosen to be i — 55°. Top: isophotes in the central 
8 pc. Bottom: cuts along (left) and perpendicular (right) 
to the P1-P2 axis. 



tral mass, self-gravity is sufficient to counteract the differ- 
ential precession of the disk. An external perturbation can 
excite this mode, and it is then long-lasting, over 100 Myr, 
or 3000 rotation periods. The interval between two such 
external perturbations (either passage of a globular clus- 
ter, or a molecular cloud) in M 31 is of the same order 
of magnitude, so that the external perturbations are an 
attractive mechanism. Each episode of to = 1 waves will 
heat the disk somewhat, but the instability is not too sen- 
sitive to the initial radial velocity dispersion. Over several 
10 8 yr periods, the nuclear disk could be replenished by 
fresh gas from the large-scale M 31 disk and subsequent 
star formation. 

6. 5. Comparison with linear WKB predictions 

In the potential of a point mass, orbits are exactly keple- 
rian, the precession rate SI — k of eccentric orbits is zero. 
The presence of a small disk of mass Md, lighter than the 
central point mass Mbh, makes this precession rate nega- 
tive, with amplitude varying as J^- . If self-gravity has 
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Fig. 21. Comparison between the STIS/HST velocity pro- 
file (crosses with error bars) and the corresponding profiles 
from an N-body simulation of an m = 1 mode where the 
disk was 40% of the central mass (with and without the 
addition of particles in the central hole; dotted and thick 
solid lines respectively). 

a large enough role, and in particular, if the disk is cold 
enough and its Jeans length smaller than the disk radius, 
m = 1 density waves can propagate; their dispersion rela- 
tion has b een st udied in t he W KB approximation (Lee & 
Goodman 1999 , Tremaine 2001 ) . In the linear approxima- 
tion, the pattern speed for the wavelength A = 2ir/k is in 
first approximation for « 1: 



TiGT, d \k\ _fcV 



n 



n 2 



where is the surface density of the disk, and T the usual 
reduction factor that takes into account the velocity dis- 
persion c of the stellar disk, and its corresponding velocity 
distribution (e.g. Tremaine |200l| ). The pattern speed then 

for a sufficiently cold disk, 



remains of the order of 



and is much smaller than the orbital frequency. These slow 

2 

waves exist whenever the thin-disk Jeans length Aj — ^jg- 
is lower than 4r, while the Toomre parameter Q is less rel- 
evant (Lee & Goodman 1999). 

Although the waves observed in the N-body simula- 
tions were non-linear, unwound and far from the WKB 
approximation, it is of interest to compare the values 
of observed pattern speed with these predictions. Since 
our disk has significant velocity dispersion, the expected 
modes are of the "pressure" nature, including the effects 
of dispersion and gravity, corresponding to the p-modes 



of Tremaine (2001). We therefore expect positive pattern 



speeds. If the effect of the background bulge can be ig- 
nored, our model simplifies into a Kuzmin disk with a 
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Fig. 22. Eccentricity e(a) of particles as a function of their 
major axis a (top); e(a) is counted positive when the 
apocentre is in the region of the PI accumulation, neg- 
ative otherwise. Polar angle (from the BH) 9 apo of the 
apocentre of orbits (middle), and 8 per i of their pericen- 
tre (bottom), at epoch 28.8 Myr. The position of the 
maximum density (corresponding to the PI component of 
M 31) is indicated. 



central point mass. It is more difficult to determine the 
amount of equivalent softening, mimicking the effects of 
velocity dispersion. A reasonable approximation is that 
the Mach number A4(r) 



VLrjc ~ 3 (see Fig. |16|). The 
reduction factor J-(kc/£l) can be approximated as a de- 
creasing exponential of krj3, where (3 ~ 1/A4 ~ 0.3. Fig. 5 



from Tremaine (2001) then predicts for the highest pattern 
speed (with the least nodes) a value of fl p ~ 2.9 km/s/pc, 
for a disk of 24% the mass of the SBH (of 7 10 7 M ), al- 
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most coincident with what we observe in the simulations 
in Fig. [I?] . This is a remarkable agreement, given all the 
approximations. Note that the Jeans length Aj maximises 
at 1.5r over the disk. 



6. 6. Excitation of the waves 

The waves can be maintained during a long time-scale, 
once excited, but triggering mechanisms remain to be 
found. The displacement of the centre of mass (and cor- 
responding unstable oscillations) is not a good candidate, 
since the frequencies are much higher. 

Tremaine (1995) has proposed that the to = 1 wave 
amplifies through dynamical friction on the bulge, if its 
pattern speed is sufficiently positive. This amplification 
results from the fact that the friction decreases the energy 
less than the angular momentum. The orbits with less and 
less angular momentum are more and more eccentric, and 
the m = 1 mode develops. To quantify numerically the dy- 
namical friction of the bulge would require a prohibitive 
number of particles in a live bulge. We have attempted 
to compute the amplitude of the effect through a semi- 
analytic method, as follows. During a simulation where 
only the nuclear disk is represented by particles as above, 
the amplitude of a possible m — 1 component, with its 
phase, are computed through Fourier transform analysis 
at each radius, and each time step. The corresponding pat- 
tern speed ttp is derived. The eccentricity of each orbit, 
averaged over 100 time-steps, is also derived. To the parti- 
cles participating in the to — 1 (with the highest eccentric- 
ities), is applied a dynamical friction force, as modelised by 
Weinberg ( 1985 ) for bars; this gives a torque proportional 
to the pattern speed f2 p , since the bulge stars are assumed 
without any significant rotation. The cumulative effects of 
the torques over several thousands dynamical times were 
found negligible. Note that in spiral galaxies like M 31, the 
bulge is slightly rotating, which can counteract the effect 
on a slow pattern rotation. 

Finally, an external perturbation is more likely to trig- 
ger the m = 1 perturbation: e.g. interstellar gas clouds 
are continuously infalling onto the nucleus. Also, the slow 
modes pattern speed discussed here are quite large with 
respect to the external disk frequencies, and resonances 
are likely to occur. 



suprisingly homogeneous throughout the nucleus after 
the bulge contribution has been subtracted, confirming 
previous results, as e.g. emphasized in KB99. However, 
contrarily to the analysis of Lauer et al. (1998), we 



7. Conclusion 

In this paper, we have first tried to clarify a number of 
issues concerning the morphology and colours of the nu- 
cleus of M 31. The photometry in the nuclear region of 
M 31 can thus be decomposed in 3 distinct components: 

— The bulge: a hot slowly rotating system with central 
colour values of {V - I) B = 1.26, (U-V) B = 1-98 

— The nucleus: exhibiting two peaks in the I band, the 
so-called PI and P2, both being offcentred with re- 
spect to the centre of the bulge isophotes. Colours are 



found that the nucleus is redder than the bulge in both 
V-I and U—V with (V-I) N = 1-30, (U-V) N = 2.36. 

— A central blue excess, the so-called UV peak, mostly 
apparent in the F300W/WFPC2 image, where it is 
just resolved. This comp onent was first detected and 
discussed by King et al. ( 1995| ). It has an axis ratio of 
about 0.7, a major-axis FWHM of ~ 0'.'21 and a PA of 
~ 62° ±8. 

We have then presented new OASIS and STIS data 
which complete our knowledge of the complex kinematics 
within the central 2" of M 31. Taking the centre of the UV 
peak in the F300W WFPC2 image as our reference zero 
point, we have derived the respective offsets required to 
reconcile the FOC data of Sta+99, the SIS data of KB99, 
and the newly presented STIS and OASIS kinematics. The 
main results are summarized here: 

— The kinematic axis (PA= 56.4°) is nearly coincident 
with the major-axis of the nucleus but not with the 
P1-P2 axis (PA= 42°). 

— The strong asymmetry of the velocity profile along the 
P1-P2 axis observed by Sta+99 is confirmed. On the 
PI side, the maximum velocity in the STIS kinematics 
(201 ±5 kms -1 ) is reached at 0'.'63 (the turnover being 
at a radius of 0742 with V — 197± 5 kms -1 ), while on 
the other side the minimum velocity(— 292±21 kms -1 ) 
is at -0'.'235. 

— The velocity dispersion peak is offset from the UV 
peak on the anti-Pi side. It reaches 321 ± 33 kms -1 
at — 0'.'235 in STIS measurements. 

— The velocity dispersion peak is consistent with being 
the result of pure velocity broadening, and is not con- 
sistent with a hot system located at the location of the 
UV peak. High velocities seem to be associated with 
P2 which is observed to be offset from the UV peak on 
the anti-Pi side. 

We finally performed new N-body simulations of to = 1 
modes in a thin disk including a SBH in its centre. We 
have found that such to = 1 modes can be maintained 
from more than a thousand dynamical times, when the 
disk mass represents between 20 and 40% of the central 
mass, allowing self-gravity to counteract the differential 
precession of the disk. The pattern speed of this m = 1 
mode is small, of the order of 3 km/s/pc in a case designed 
to resemble M 31's nucleus. For the present paper, we have 
only presented a model where the m = 1 mode was excited 
by launching particles along elliptical orbits. This model 
reproduces reasonably well the observed asymmetry in the 
surface brightness distribution and mean stellar velocity, 
although we had to remove some particles near the SBH 
to avoid numerical artefacts. A detailed comparison with 
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the observed stellar velocity dispersion requires a realistic 
treatment of the particles near the SBH. The existence 
of an m — 1 mode in the nucleus of M 31 is consistent 
with the presence of a SBH with a mass in the range ~ 
3.5- 8.5 10 7 M Q . 

The model proposed here relies on a relatively thin and 
cold nuclear stellar disk. How can this disk be formed and 
maintained? A likely possibility is that the central region 
is subject to accretion and infall of gas from the nearby 
disk of M 31, which is gas rich. Already, within the few 10- 
100 pc dust lan es, an d CO molecular clouds are observed 
(Melchior et al. 2000 ). The infall of gas clouds could both 
trigger the m = 1 instability, and also progressively reform 
the cold stellar disk. 
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Appendix A: Ionised gas in the centre of M 31 

We have reexamined data obtained with the TIGER spec- 
trograph (CFHT) in the spectral domain around the 
[Nn]A6548/Ha A6563 emission lines, already discussed in 
Bacon et al. (1994), to more systematically search for the 
possible presence of ionised gas. The detection of a po- 
tential Ha emission line is difficult due to the underlying 
deep stellar Ha absorption feature. We have thus applied a 
new algorithm (Emsellem et al., in preparation) allowing 
an accurate subtraction of the stellar continuum, which 
exploits a complete library of stellar and galaxy spectra 
with 0.5 A sampling. The typical detection limit (3cr) is^] 
then 3.6 10 -19 W.m -2 . 

We are now successfully detecting an spectrally unre- 
solved source of Ha A6563 line emission, about 2 / / 43 ( C/39 
North, 2"4 East) from the centre of M 31 (Figs. |Aj] and 
|A.2| ) . This emission line was previously undetected mainly 
because it is well sunk into the corresponding absorption 
feature, The emission intensity peaks at more than 5 times 
above the noise with a redshift of — 45 ± 10 kms -1 (helio- 



centric, see Fig. AT): this is to be compared with the sys- 
temic velocity of —308 kms -1 , as derived from the OASIS 
spectra, which gives a relative velocity of 263 kms -1 . This 
emission line is observed on several adjacent spectra, its 
distribution being consistent with a spatially unresolved 
system at the resolution of the TIGER data (FWHM of 
~ l'/l). This is not an instrumental artefact, as neigh- 
bouring spectra come from non adjacent regions of the 
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Fig. A.l. Top: original spectrum of M 31, summed over a 
region where Ha emission is found (5 spectra). The dotted 
line shows the fit performed using a library of high resolu- 
tion spectra broadened and shifted according to the local 
kinematics. Bottom: resulting spectrum after subtraction 
of the template fit. The unresolved Ha A6563 line clearly 
stands out, more than 5 sigma above the noise. 



CCD and exhibit the same spectral line. This is further 
supported by the detection of a point-like source in the 



continuum subtracted F658N/WFPC2 image (see Fig. A. 2), 
as well as by recently published two-dimensional spec- 
troscopy (del Burgo, Mediavilla and Arribas 2000] ; source 
"A"). Note that the central sources A, B, C and D dc- 



Note that the detection limit quoted in Bacon et al. ( 1994 ) 
was wrong due to an error in the flux units. 



tected by del Burgo et al. ( [2000D all appear as point-like 
in the F658N/WFPC2 image, their fluxes being consistent 
with single planetary nebulae, in contradiction with the 
claims made by these authors. 

This feature can be simply explained by a single plan- 
etary nebula. Indeed the integrated flux of this point- 
like source is 12.5 ± 3 10 -19 W.m -2 in the WFPC2/F658N 
continuum subtracted image (we derive a lower limit of 
10 10 -19 W.m -2 from the TIGER datacube), compatible 
with a planetary in M 31. The velocity of the detected 
Ha emission line is consistent with the expected stellar 
velocities of ~ 250 kms -1 at the edge of the nucleus on 
the PI side, but is not consistent with the low velocities 
(< 20 kms -1 ) of bulge stars. This strongly suggests that 
the object associated with this emission line belongs to the 
nucleus itself. If we now assume that this emitting source 
is in the plane of the nucleus, it is then at ~ 24° from the 
line of nodes, and at 2'.'81 (or 10.5 pc) from the UV peak 
(for an assumed inclination of i = 55°). 
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